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1 Introduction 

So far, there exist many explosion models proposed by theorists for different types of supernovae, but still 
there are no definite criteria to decide which of the models are realized in nature. Only a few parameters, 
such as kinetic energy and total ^^Ni production, can be derived directly from the modelling of the 
explosion and compared with the observational values. The subsequent evolution of the exploded star 
gives us much more possibilities to compare models and to decide which one fits observations better. 

The first possibility is modelling of 7-ray luminosity of supernovae. The comparison with the obser- 
vational values allows to define the total mass of radioactive isotopes, to judge on the composition of the 
outer layers where 7-rays are absorbed and thermalized, and also to check the approximations of 7-ray 
opacity. During the first months after an explosion one can examine a theoretical model by calculating 
bolometric and monochromatic light curves and spectra, and comparing them with observations. Later 
on, gas in the ejecta cools down and becomes almost unobservable. The next opportunity to analyze the 
ejecta is on the stage of a young supernova remnant (SNR), when noticeable amount of circumstellar gas 
is swept up. Then a reverse shock forms, goes inwards the ejecta and illuminates them once again. 

The successful theoretical supernova explosion models should be able to explain any feature of the 
emission from supernovae at any evolutionary stage. 

Different combinations of codes we have in our group allows us to compute the evolution of models, 
and therefore test them, on several stages. For calculations of broad-band optical and bolometric light 
curves we use the multi-frequency radiation hydro code STELLA. At the young remnant stage, while gas 
is transparent, we calculate the evolution and X-ray emission by combination of hydro part of STELLA 
with the non-equilibrium ionization code which is based on the original algorithm by Peter Lundqvist. 



2 SN la models 

We choose five models of Type la supernovae which are shown in Table 1. W7 [1] and DD4 [2] are more or 
less similar classical ID Chandrasekhar mass models. LA4 [3] and WD065 [4] are sub-Chandrasekhar mass 
ones. (See original works for details.) The last model was computed a couple of years ago in Garching, in 
the Max-Planck-Institute group [5]. It is originally 3D. We average it over An for our calculation, since 
we use only ID codes at the moment. The main feature of the model is that it is very well mixed, while 
energetics and the amount of Ni^^ are lower than in other Chandrasekhar mass models. In our calculations 
we compare two versions of the model with the same hydro part, but different element distributions. 



3 Results 
3.1 SNe la 

Now there appears more and more possibilities to obtain 7-ray light curves using modern 7-ray space 
observatories, like COMPTEL [6]. Here we just compare what one would observe after the explosion of 
the models we have discussed. 
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Table 1. Parameters of SN la models 



Model DD4 W7 LA4 WD065 MR 



MwD^ 1.3861 1.3775 0.8678 0.6500 1.4 

MssNi'' 0.63 0.60 0.47 0.05 0.42 

Esi^ 1.23 1.20 1.15 0.56 0.46 

''in Mq 

''in lO''^ ergs s"^ 
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Fig. 1. 7-ray production by Ni and Co decays (dashes), 7-ray power deposited into heating (dots), and total 
7-ray luminosity of SN (solid) versus time in days for different explosion models: W7 {top panel), MR model 
averaged over 47r (bottom left), and MR model averaged over an opening angle of 14° with the same mass of 5^Ni, 
which is situated mostly near the center (bottom right). 
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In fig. 1 we compare the light curves in 7-rays produced by different models. We have chosen W7 as 
a representative case of Chandrasckhar-mass models. The light curves in 7-rays arc more or less similar 
for all of them. The luminosity of Sub-Chandrasekhar-mass MR model is lower due to smaller amount of 
^^Ni, and it rises slower since it is less energetic, so 7-photons are locked inside the ejecta for longer time. 
One can sec even smaller mimbcr of 7-photons during the first weeks after the explosion in the model 
similar to MR but less mixed, with ^^Ni residing mostly near the center. The maximum luminosity in 7 
is still the same as in original MR, since we preserve the total amount of ^^Ni. 

Therefore, the light curve in 7 rays during the first months can be divided into two epochs, that 
represent different physical parameters of the explosion: the first 40-60 days and the following evolution. 
From the observations during the first period one can judge on the combination of explosion energy and 
the distribution of radioactive stuff over the ejecta, while the second period tells us mostly about total 
production of radioactive elements at the explosion. Since currently the observations can only provide us 
with fluxes integrated over several days (or even weeks), it seems expedient to make two different sets 
of observations: during the first two months after the explosion, and after that, in order to distinguish 
between different parameters of explosion, and to test explosion models. 

We have described our modelling of bolometric and UBVRI light curves and calculations of X-ray 
emission of young SNRs in the previous work. For details one can see e.g. [7, 8, 9]. We just mention here 
that a bit less energetic and more mixed models, like MR, seem to us more preferable both for SN light 
curve and for SNR X-ray emission. At the latter point, our results seem to contradict the conclusions 
by [10], that iron lines from well mixed models at the SNR stage arc too strong, while less mixed models 
give a better agreement with observations of Tycho. Still we trust in our results more, since we take into 
account the ionization energy in the equation of state, which is comparable to the thermal energy. Electron 
thermal conductivity is also included into our calculations, and it smoothes very much the temperature 
profile between forward and reverse shocks. Both these effects are neglected in the code by [10], though 
they are able to change the emitted spectrum appreciably. The code by [10] takes into account possible 
difference between electron and ion temperature using the standard Coulomb collision equation, while we 
treat this effect parametrically. This can also lead to differences in our results. The work on improving 
physics and making our code more self- consistent is in progress. 

3.2 SN 1993J 

In order to pay attention to the SN entitling this conference, we have recalculated its light curve. We use 
the same model 13C from [11] which was investigated a few years ago in [12]. It is one of the best models 
for SN 1993J, and we just wanted to check if the new version of our code STELLA, with improved and 
renovated physics, still produces the light curve that fits observations well. 

The results of previous and current calculations are compared in fig. 2. The main improvement in the 
code is a new approach for expansion opacities [13]. They become much more complicated, and, most 
probably, the bumps on the new light curves are the results of this improvement. But sometimes (for 
instance, in the U band) these bumps seem to fit the observations even better than it was in the old 
version. R and I bands are calculated with STELLA for the first time, and they look perfect. 

There is a discussion in the literature on the estimates of the reddening to SN 1993 J ([14] and references 
therein). In the plots in fig. 2 we assume E{B — V)= 0.08. With larger value, the observational curves 
become higher than the modelled ones. So the new version of STELLA confirms that the model 13C 
corresponds to the explosion of SN 1993J very well, but in the case of low reddening. 
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Fig. 2. Bolometric and UBVRI light curves for SN 1993J. Top left pane? shows the best model 13C8 from [12] {solid 
line - calculations with STELLA, dashes - with EDDINGTON), top right and bottom panels: our recalculation 
of the same model with a new version of STELLA. Observational data are taken from [14] and shown by crosses. 
The value of reddening is taken to be E{B — V) — 0.08. 



7. S. Blinnikov, E. Sorokina: Type la Supernova models: latest developments. In: Hunting the cosmological 
parameters, ed by D. Barbossa et al (KLUWER 2003) (astro-ph/0212530) 

8. S. Blinnikov, E. Sorokina: A&A 356, L30 (2000) 

9. D.L Kosenko, E.LSorokina, S.L Blinnikov, P. Lundqvist: X-ray emission of young SN la remnants as a probe 
for an explosion model. In: Proc. of 34th Scientific Assembly of COSPAR, (ASR 2003) (astro-ph/0212188) 

10. C. Badenes, E. Bravo, K.J. Borkowski, I. Dominguez: ApJ 593 358 (2003) 

11. S.E. Woosley, R.G. Eastman, T.A. Weaver, RA. Pinto: ApJ 429, 300 (1994) 

12. S.I. Blinnikov, R. Eastman, O.S. Bartunov et al: ApJ 496, 454 (1998) 

13. E.I. Sorokina, S.I. Blinnikov: Energy exchange inside SN ejecta and light curves of SNe la. In: Nuclear 
Astrophysics, MPA/P13, ed by W. Hillebrandt, E. Miiller (Garching 2002) pp 57-62 

14. M.W. Richmond, R.R. Treffers, A.V. Filippenko et al: Astron. J. 107, 1022 (1994) 



